Spin-flip scattering in time-dependent transport through a quantum dot: 
Enhanced spin-current and inverse tunneling magnetoresistance 
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We study the effects of spin-flip scatterings on the time-dependent transport properties through 
a magnetic quantum dot attached to normal and ferromagnetic leads. The transient spin-dynamics 
as well as the steady-state tunneling magnetoresistance (TMR) of the system are investigated. The 
absence of a definite spin quantization axis requires the time-propagation of two-component spinors. 
We present numerical results in which the electrodes are treated both as one- dimensional tight- 
binding wires and in the wide-band limit approximation. In the latter case we derive a transparent 
analytic formula for the spin-resolved current, and transient oscillations damped over different time- 
scales are identified. We also find a novel regime for the TMR inversion. For any given strength of 
the spin-flip coupling the TMR becomes negative provided the ferromagnetic polarization is larger 
than some critical value. Finally we show how the full knowledge of the transient response allows 
for enhancing the spin-current by properly tuning the period of a pulsed bias. 

PACS numbers: 73.63.-b,72.25.Rb,85.35.-p 



I. INTRODUCTION 

Single-electron spin in nanoscale systems is a promis- 
ing building block for both processors and memory stor- 
age devices in future spin-logic applications; 1 - In partic- 
ular great attention has been given to quantum dots 
(QDs), which are among the best candidates to imple- 
ment quantum bit gates for quantum computation! 2 - In 
these systems the spin coherence time can be orders of 
magnitude longer than the charge coherence timeSf 3 -^^ 
a feature which allows one to perform a large number of 
operations before the spin-coherence is lost. Another ad- 
vantage of using QDs is the possibility of manipulating 
their electronic spectrum with, e.g., external magnetic 
fields and gate voltages, and fine-tuning the characteris- 
tic time-scales of the system. 

For practical applications, it is important to achieve 
full control of the ultrafast dynamics of QD systems af- 
ter the sudden switch-on of an external perturbation. 
To this end the study of the time evolution of spin- 
polarized currents and the characterization of the spin- 
decoherence time is crucial. The theoretical study of the 
transient regime is also relevant in the light of recent 
progresses in the time-resolution of dynamical responses 
at the nanoscale. Novel techniques like transient current 
spectroscopy^ and time-resolved Faraday rotation 8 -^ al- 
low us to follow the microscopic dynamics at the sub- 
picosecond time-scale. These advances may open com- 
pletely new scenarios with respect to those at the steady- 
state. For instance, transient coherent quantum beats of 
the spin dynamics in semiconductor QD's have been ob- 
served after a circularly polarized optical excitation,-^ 
and theoretically addressed by Souza.— We also mention 
a recent experiment on split-gate quantum point contacts 



in which the measured quantum capacitance in the tran- 
sient regime was six orders of magnitude larger than in 
the steady-state.— 

Another attracting feature of magnetic QDs is their 
use as spin-devices in magnetic tunnel junctions (MTJs). 
Different orientations of the polarization of ferromag- 
netic leads results in spin-dependent tunneling rates and 
hence in a nonvanishing tunneling magnetoresistance 
(TMR). In recent experiment a 12 i 13 ' 14 the inverse TMR 
effect (TMR < 0) has been observed and various models 
have been proposed to address such phenomenon ! 15 ' 16 ! 17 

In this paper we study the time-dependent transport 
through a single level quantum dot connected to nor- 
mal and ferromagnetic leads. In order to get a sensi- 
ble transient regime, we adopt the so-called partition- 
free approach, in which the electrode-QD-electrode sys- 
tem is assumed to be in equilibrium before the external 
bias is switched on . 18 ' 19 Spin-symmetry breaking terms 
(like spin-flip scatterings and spin-dependent dot-leads 
hoppings) are included— and this requires the time- 
propagation of a genuine two-component spinor. 

Explicit calculations are performed both in the case 
of non-interacting one-dimensional (ID) leads of fi- 
nite length and in the wide-band-limit approximation 
(WBLA). In the WBLA we derive a closed analytic for- 
mula for the exact spin-resolved time-dependent cur- 
rent, which can be expressed as the sum of a steady- 
state and a transient contribution. The latter consists 
of a term which only contains transitions between the 
QD levels and the electrochemical potentials {resonant- 
continuum), and a term which only contains intra-dot 
transitions [resonant-resonant). Remarkably these two 
terms are damped over two different time-scales, with 
the resonant-continuum transitions longer lived than the 
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resonant-resonant transitions. We further show that, go- 
ing beyond the WBLA, extra transitions occur. These 
involve the top/bottom of the ID electrode bands and 
might be relevant for narrow-band electrodes. 

For QDs connected to normal leads we study the quan- 
tum beat phenomenon in the presence of intra-dot spin- 
flip coupling of strength V s f . We show that the amplitude 
of the beats is suppressed independently of the structure 
of the leads (WBLA and ID leads). In addition we show 
how to engineer the spin-polarization of the total current 
by exploiting the full knowledge of the time-dependent 
response of the system. By applying a periodic pulse 
of proper period, we provide numerical evidence of os- 
cillating spin-polarizations with amplitude two orders of 
magnitude larger than in the DC case. 

Finally, in the case of ferromagnetic electrodes we 
study the steady-state TMR and the conditions for its 
inversion. Treating the leads in the WBLA we derive 
a simple formula for the TMR in linear response. No- 
ticeably, in the presence of spin flip interactions there is 
a critical value of the polarization of the ferromagnetic 
leads above which the TMR is negative even for sym- 
metric contacts to the left /right leads. This provides an 
alternative mechanism for the TMR inversion. The above 
scenario is qualitatively different with ID leads, e.g., the 
TMR can be negative already for V s f = 0. 

The paper is organized as follows. In Section |TT] we 
introduce the lead-QD-lcad model. In Section IIIII we 
employ the WBLA and derive an exact formula for the 
time-dependent current in the presence of spin-symmetry 
breaking terms. We present explicit results both for nor- 
mal and ferromagnetic leads, focussing on transient quan- 
tum beats (normal) and the TMR (ferromagnetic). In 
Section IIVI we address the above phenomena in the case 
of ID leads and engineer the dynamical spin responses. 
Here we also describe the numerical framework employed 
to compute the time-dependent evolution of the system. 
The summary and main conclusions are drawn in Section 

II. MODEL 

We consider the system illustrated in Fig Q] which con- 
sists of a single-level QD contacted with ID Left (L) and 
Right (R) electrodes. The latter are described by the 
tight-binding Hamiltonian 

N 

H a = (^ajaCaja-Caja 

+ VoC^Caj+ic + V *cl J+llT C aja ^ , (1) 

where N is the number of sites of lead a = L,R, the 
nearest-neighbor hopping integral is Vq and c^J a is the 
annihilation (creation) operator of an electron on site 
i of the lead a with spin a. In the case of ferromag- 
netic leads, we distinguish between two configurations, 



one with parallel (P) and the other with antiparallel (AP) 
magnetization of the two leads. In order to model these 
two different alignments we use the Stoner prescription 
according to which the on-site energies are 



lead I QD lead J? P (APj 
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FIG. 1: (Colour online) Schematic illustration of the model 
Hamiltonian. The lead L with majority spin f electrons and 
the lead R with majority spin f (J.) electrons in the P (AP) 
configurations are separated from the QD by tunnel barriers 
which are accounted for with a renormalized hopping V . The 
energy-spin-splitting is 2Ae in the leads and 2E Z in the QD. 
Spin- flip scatterings of strength V s f occur in the QD. 



e aja = (-l)^Ae P 

e a]a = (-lj'-'C-l^Ae AP (2) 

where 2Ae is the band spin splitting. In the above model 
the P configuration corresponds to a majority-spin elec- 
tron with spin f in both leads, while in the AP configura- 
tion the majority-spin electron have spin \ in lead L and 
spin J, in lead R. The case of normal leads corresponds 
to Ae = 0. The Hamiltonian of the quantum dot reads 

Hqd = £ <iod\d a + V.fC^di + K)4 d T , (3) 

a 

where annihilates (creates) an electron on the QD 
with spin a and V s t is responsible for intra-dot spin- 
flip scatterings. The on-site energy £<j CT = £<i + &E z /2, 
fj = ±1, where E z is an intra-dot energy splitting. The 
quantum dot is connected to the leads via the tunneling 
Hamiltonian 

h t = E {Vdo, a o>d\c ala , + V^ aa ,cl la ,d^ . 

ck—L,R aa' 

(4) 

with Vda,aa' the amplitude for an electron in the QD 
with spin a to hop to the first site of lead a with spin 
a'. Alternatively, one can express Ht in terms of the 
one-body eigenstates of H a labelled by (a, k, a), 

N 

ct=L,R k—1 oa' 

(5) 

with Vda.kaa' = ^ 2 / (N + 1) V dc7 , a a> sink. 

Putting all terms together the Hamiltonian of the 
whole system in equilibrium reads 

H Q =H L +H R + H QD + H T . (6) 
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In the next two Sections we study out-of-equilibrium 
properties of the system described in Eq.©, after a sud- 
den switch-on of a bias voltage U a (t) in lead a = L,R. 
We calculate the time-dependent spin-polarized current 
I aa {t) flowing between the QD and lead a. I Qcr (i) is de- 
fined as the variation of the total number of particles N aa 
of spin a in lead a, 



Iaa(t) = -JT^aa — 23? ^ Vako.do' G\ 
M k.o> 



a 1 .Oika 



(*,*), (7) 



where G^ a , aka is the dot-lead lesser Green's function of 
the contacted non-equilibrium system, and 5ft stands for 
the real part. Unless otherwise stated, the current is 
computed at the L interface and the short-hand nota- 
tion I /y — lL a is used for the spin-polarized current while 

= If — 1\ denote the total cur- 



Itot — If + h and I spi7 
rent and spin current respectively. We further specialize 
to constant biases U a (t) = U ai except in Section HV Bl 
where the bias is a periodic pulse. Therefore we define 
U = Ul, while the value of Ur is either or —U. The nu- 
merical simulations have been performed using dot-leads 
hoppings independent of a (symmetric contacts). We 
show that, despite the symmetry between L and R leads, 
the sign of the steady-state TMR can be inverted, a prop- 
erty which has been so far ascribed to strongly asymmet- 
ric tunneling barriers. 15 ! 21 Moreover we specialize Eq.Q 
to the case V d ], a ] = Vdi, a { = V, and V d -\, a i = Vdi,at = 0- 
However we note that spin-flip scatterings are accounted 
for by considering a finite intra-dot V s f which is taken to 
be real in this work. If not otherwise specified, the initial 
Fermi energy Ep of both leads is set to (half-filled elec- 
trodes). In the numerical calculations below all energies 
are measured in units of Vq (4Vq is the bandwidth of the 
leads), times are measured in units of H/Vq and currents 
in units of eVo/h with e the electron charge. 



QD hops to lead a by changing its spin from a to is, and 
hops back to the QD with final spin a' . Exploiting the 
Dyson equation the expression of (S aj/ ) fflT / is 



V / ( 



E 

k 



v d (OJ) Vakv.da' , (9) 



where g a kv = (u> — e a k v — U a + ir])^ 1 is the retarded 
Green's function of the uncontacted a lead expressed 
in term of the one-particle eigenenergies e a ku- In the 
WBLA Eq.© becomes 



(E a „H)^-j(l 



(10) 



i.e., the self-energy is independent of frequency. The 2x2 
matrices T's have the physical meaning of spin-dependent 
tunneling rates and account for spin-flip processes be- 
tween the leads and the QD. We wish to point out that 
each r matrix must be positive semi-definite for a proper 
modeling of the WBLA. Indeed, given an arbitrary two 
dimensional vector («| , V[ ) one finds 

^V*(T qi ,) V a > — 2-77^ 5(u)-£ a kv-Ug)\ ^ Vda.akyVg 
a. a 1 k a 

(11) 

The above condition ensures the damping of all transient 
effects in the calculation of local physical observables. 

Proceeding along similar lines as in Reffl9l one can 
derive an explicit expression for spin-polarized current 
I aa (t) defined in Eq.|7]), 



(12) 



III. RESULTS IN THE WIDE BAND LIMIT 
APPROXIMATION 



f(uj) being the Fermi distribution function. In the above 
equation I^ a is the steady-state polarized current which 
for a = L reads 



Here we study the model introduced in the previous 
Section within the WBLA. This approximation has been 
mostly used to study spinless electrons and a closed 
formula for the exact time-dependent current has been 
derived i ' Below we generalize these results to sys- 
tems where the spin symmetry is broken. The retarded 
Green's function projected onto the QD can be expressed 
in terms of the embedding self-energy 



£(w) 



j/=T,J. a=L,B, 



which is a 2 x 2 matrix in spin space. In Eq. JH]) (Y, av ) aa i 
accounts for virtual processes where an electron on the 



La 



2^ 



x Tr 



(/(w - U L ) - f(cu - U R )) 



1 



1 



uj - h qd + f r u- Hqd - f r 



(13) 



with r Q = J2 a Tao- an d r — J2 a ^a- The steady current 
at the right interface I Ra is obtained by exchanging L 



(8) R in the r.h.s. of Eg-lp]), We observe that I Ra ± -I[ 



since the spin current is not conserved in the presence of 
spin-flip interactions. The second term on the r.h.s. of 
Eq. p2|) describes the transient behavior and is expressed 
in terms of the quantity 
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T a (uj,t) = 



5> 



e i(nj+U fJ -H QD + ± r)t 



id, 



1 



Q,/3 



Up - Hq D + |] 
e i(ffQD + |r)t 



(w - Hqd + |r)( w + - #qd + |r) 

(lu - H QD + ir)(o; + - Hqd + \ f) 1 " (w - Hqd - ±r)(w + - Hqd 



(14) 



Few remarks about Eq. (|14p are in order. In linear re- 
sponse theory only the contribution of the first line 
remains. Such contribution is responsible for tran- 
sient oscillations which are exponentially damped over a 



timescale | tj | and have frequencies ui c 



\E F + U a 



where hj + ir^ , j = 1,2, are the two eigenvalues of 
Hqd + iT/2. These oscillations originate from virtual 
transitions between the resonant levels of the QD and the 
Fermi level of the biased continua (resonant- continuum) . 
No information about the intra-dot transitions (resonant- 
resonant) is here contained. Resonant-resonant transi- 
tions are instead described by the contribution of the 
second line of Eq.(|14p which yields oscillations of fre- 
quency u>x t 2 — \h\ — hq\ damped as exp(— t/ri 2) with 
= r i| _1 + l T 2| _1 - Finally, it is straightforward to 
verify that for T a and Hqd both diagonal in spin space, 
Eq. (fT2f decouples into two identical formulae (but with 
different parameters), which exactly reproduce the well- 
known spinless result . 19 i 22 

Below we consider diagonal matrices T avi in agreement 
with the discussion at the end of Section [TTJ Thus, all 
spin- flip scatterings occur in the quantum dot. The ferro- 
magnetic nature of the leads is accounted for by modeling 
the r matrices as 2 -^ 



a = L, R 
a = L.R 



AP 





TT 


= r 0Q (i 


+ p) 




U 


- r a(i 


-p) 




rr 


= r oz (i 


+p) 


tli ] 


u 


= r 0L (i 


-p) 




'it 


= r ofl (i 


-p) 




' u 


= r «(i 


+ p) 



(15) 



where all the remaining matrix elements are vanishing 
and < p < 1 is proportional to the polarization of 
the leads. Moreover we assume that To a = Tq does not 
depend on a, yielding a Left /Right symmetry in absence 
of ferromagnetism. 



A. Normal case: quantum beats 

In this Section we consider a quantum dot with an 
intra-dot energy splitting E z between the f and j levels, 
and an intra-dot spin- flip energy V s t. The QD is coupled 



to L and R normal electrodes, p — 0. The spin split- 
ting E z leads to two different transient frequencies and 
produces coherent quantum beats in both the total and 
spin currents^ The effect of spin-flip scatterings on the 
coherent oscillations is studied. 

To make contact with Reffiol we consider the same 
parameters, i.e. Ul = 0, Ur = 200, Eda — &(t)Un/2 + 
aE z /2 with E z = 10, Tq = 1 and inverse tempera- 
ture (3 — 100. However the intra-dot spin-flip coupling 
V s f is, in our case, non zero. In Fig|5] we show the 
time-dependent current Itot(t) and the discrete Fourier 
transform^ I tot (uj), of Itotif) — hotipo) for different val- 
ues of V s f — 0, 10, 20. The spin current I sp i n (t) is dis- 
played in Fig[3J The small tunneling rate T leads to a 
very long transient regime, a property which allows us to 
observe well defined structures in the Fourier spectrum of 
the current. At V s f = the frequencies of both I^(t) and 
Il(t) are w a \ = uj a 2 = \Ur/2±E z /2\ — 95, 105, in agree- 
ment with the results of Reffiol. This is confirmed by the 
Fourier analysis of I to t which is displayed in panel b) of 
Figl^l As V s f is increased the frequencies of the transient 



oscillations change according to \Ur/2 ± ^JE^/A + V^\, 

see panels d) and f) of FigJ^l and the amplitude of the 
quantum beats in I to t and I sp i n is suppressed. Such sup- 
pression is due to the fact that the Hamiltonian of the 
whole system is spin-diagonal if we choose the quantiza- 
tion axis along 



£ = I cos 9 + x sin . 



(16) 



with cos 6 = E z j 'JE% + This follows from the in- 

variance of the normal electrode Hamiltonians (p = 0) 
and of the tunneling Hamiltonian (T aa = Fol2x2) under 
rotations in spin-space. Therefore, the time-dependent 
spin current measured along z is smaller the larger V s f 
is. Of course the spin current measured along the quan- 
tization axis £ is not suppressed by changing V s f. 

Despite the quantum beats phenomenon and its sup- 
pression with increasing V s f are well captured within the 
WBLA, only a sub-set of transient frequencies can be 
predicted in this approximation. To describe transitions 
between the top/bottom of the electrode bands and the 
resonant levels requires a more realistic treatment of the 
leads Hamiltonian, see Section [TVl 
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FIG. 2: I to t(t) and I tot (u>) for V sf = [panels a)-b)], V s} = 10 
[panels c)-d)], V s f = 20 [panels e)-f)]. The discrete Fourier 
transform is calculated using 50 equidistant points of Itot(t) — 
Itot(oo) with t in the range (0.5,1)- The rest of the parameters 
are U L =0,U R = 200, e drT = 6(t)U R /2+aE z /2 with E z = 10, 
To = 1 and inverse temperature /3 = 100. 
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FIG. 3: (Colour online) I sp in(t) for V sf = 0, 10, 20 and the 
same parameters as in FigO 



B. Ferromagnetic case: TMR 

The finite spin-polarization p of the electrodes breaks 
the invariance of Hi and Hr under rotations in spin- 
space. Thus the one-particle states become true spinors 
for nonvanishing V s f. Using Eq.(fl3| we calculate the 
steady-state total currents both in the P and AP config- 
urations, which we denote as Ip and Iap, and study the 
TMR, 



TMR: 



Ip — Iai 
Iap 



(17) 



For V s f = the TMR is always positive and its sign 
can be inverted only during the transient regime at 
sufficiently low temperatures^ The positiveness of the 
steady state TMR is due to the fact that T i = T r. It 
is possible to show that for T i 3> Toi? (or r L <C Tqr) 



and Eda — 0, the TMR w — p 2 in linear response theory 
yielding the so-called resonant inversion of the TMR. 15 
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FIG. 4: (Colour online) Contour plot of the TMR at the 
steady-state in units of 10~ 2 as a function of U and p for 
different values of V 3 f == 0.3 (top-left), 0.33 (top-right), 0.4 
(bottom- left), 0.47 (bottom right). The boundary TMR = 
is displayed with a black line. The remaining parameters are 
r = 0.5, e ia = 0, P = 100. 

In FigH] we display the contour plot of the TMR in 
the parameter space spanned by the bias voltage U = 
Ul = —Ur and the polarization p for different values of 
the spin-flip energy V s f and for E z = 0, Tq = 0.5 and 
(3 = 100. We observe that by increasing V s f, a region 
of negative TMR appears and becomes wider the larger 
V s f is. The region in which the TMR is appreciably dif- 
ferent from zero (|TMR| > 0.05 ) is confined to the high 
magnetization regime, i.e. p > 0.4. The minimum of the 
TMR is reached for the largest value of V s f at small U 
and large p, see panel d) of FigfJ] For small biases the 
sign inversion of the TMR can be understood by calcu- 
lating the currents Ip and Iap in linear response. By 
expanding Eq. (|13p to first order in U and using the T 
matrices of Eq. (Ti"5")) one can show that to leading order 
in the bias 



TMR = p 2 



(1 



p 2 n-(i- P 2 ) 



2, T 2 V 2 



[V 2 f + (1 -P 2 )T 2 }[(1 +P 2 )V 2 f + (1-P 2 )T 2 } ■ 

For any given V s f and To the denominator in the above 
expression is positive while the numerator can change 
sign. The regions of posite/negative TMR are displayed 
in the left side of FigLU For < 2V 2 f the TMR is 
negative independent of the polarization p. This follows 
from the fact that the tunneling time is larger than the 
time for an electron to flip its spin, thus favoring the AP 
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alignment. Along the boundary Tq = 2V£t the TMR is 
zero only for p = and negative otherwise. In the region 
Tq > 2Vgf there exists a critical value of the polarization, 

pc = Ji- r4 2 pL 2 , (19) 

such that the TMR is positive for p > p c and negative for 
p < p c . It is worth to emphasize that at the boundary 
V s f = the TMR is positive for any finite value of p and 
zero for p — 0. This analysis clearly shows that the sign 
of the TMR can be reversed without resorting to asym- 
metric couplings Tl i= Tr, provided spin- flip processes 
are included. The right side of Fig [5] displays the TMR 
as a function of the polarization for different values of 
the ratio V^/r 2 , < 1/2 (region of TMR inversion). 




FIG. 5: (Colour online) Left panel: in the region Tq < 2Vgt 
the TMR is always negative while in the region rj] > 2V s 2 f the 
TMR changes sign depending on the value of the polarization 
p. Right panel: TMR as a function of p for To = 1 and 
different values of V s f. 

We also have investigated the temperature dependence 
of the TMR. This is relevant in the light of practical ap- 
plications, where a large TMR is desirable at room tem- 
perature. The increase of temperature tends to suppress 
the negative values of the TMR while the positive val- 
ues remain almost unaffected. When j3 f» 1 the region of 
the TMR inversion disappears. The positiveness of the 
TMR at high temperatures has been already observed in 
ReflH. 

In the next Section we provide an exact treatment of 
the ID leads and illustrate differences and similarities 
with the results obtained so far. 



IV. RESULTS FOR ONE-DIMENSIONAL LEADS 

The numerical results contained in this Section are ob- 
tained by computing the exact time-evolution of the sys- 
tem in Eq.(|5J) with a finite number N of sites in both L 
and R leads . 25 ' 26 ! 27 Let us define the biased Hamiltonian 
at positive times as 

H bms (t)=H + H(t) (20) 



with 

N 

H(t) = 6(t) J2 E u ^yl^ a3 a . (21) 

Both Ho and Hi,i as have dimension 2(2N + 1), where the 
factor of 2 accounts for the spin. 

We use the partition-free approach 18 and specialize to 
a sudden switching of a constant bias. Accordingly, we 
first calculate the equilibrium configuration by numer- 
ically diagonalizing H$, and then we evolve the lesser 
Green's function as 

G< (t, t') = ie- iH »^ t f(H )e iHb ^ t ' , (22) 

with / the Fermi distribution function. The spin- 
polarized current flowing across the bond m — n is then 
calculated from 

■manner' 

where [. . . ] mo - ;rlo -< denotes the matrix element associ- 
ated to site m with spin a and site n with spin a', 
while 3 stands for the imaginary part. The above ap- 
proach allows us to reproduce the time evolution of the 
infinite-leads system provided one evolves up to a time 
T m ax ~ 2N/v, where v is the maximum velocity for an 
electron with energy within the bias window. Indeed for 
t > T max high-velocity electrons have time to propagate 
till the far boundary of the leads and back yielding un- 
desired finite-size effects in the calculated current, see 
FigJS] For this reason we set N much larger than the 
time at which the steady state (or stationary oscillatory 
state in the case of AC bias) is reached. We tested this 
method by comparing our numerical results against the 
ones obtained in Refs.24,28 where the leads are virtu- 
ally infinite, and an excellent agreement was found for 
t < T max . Moreover, the value of the current at the 
steady-state agrees with the Landauer formula with high 
numerical accuracy. 

Below we study the spin-polarized current at the left 
interface I a — I\h,d,a flowing between the first site of 
the lead L (m = XL) and the QD (n = d) with spin 
a. We stress that our approach is not limited to the 
calculation of the current through a specific bond, as we 
have access to the full lesser Green's function. Sensible 
electron densities and currents in the vicinity of the QD 
can be extracted and their calculation requires the same 
computational effort. 

A. Normal case: quantum beats 

We study the quantum beats phenomenon for normal 
ID leads (Ae = 0) and intra-dot spin-splitting E z for dif- 
ferent values of the spin relaxation energy V s f . The com- 
parison between the results obtained within the WBLA 
and with ID tight-binding leads will turn out to be very 
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FIG. 6: (Colour online) if(t) at zero temperature for N = 
700, U L = 0.5, Ur = 0, E F = (T max » 700), and £ F = 
— 1.7 (T maa; 870). The rest of the parameters are e da — 0, 
V = — 1, V 3 f = 0. The current reaches a well-defined steady 
state before the occurrence of finite-size effect (t > T ma x)- 



useful to elucidate advantages and limitations of the for- 
mer approach. 

In Fig[7|we show the time-dependent currents Itot{t) 
and its discrete Fourier transform Itot{&>) at zero temper- 
ature using the parameters Ul = 0.7, Un = 0, = 0.03, 
s d a = 0{t)U L /2 + <tE z /2, with E z = 0.04, for different 
values of V s f — 0, 0.05, 0.1. The spin-current for the 
same parameters is displayed in Fig[8] The value of the 
dot-lead hopping V corresponds to an effective (energy- 
dependent) tunneling rate T of the order of 10~ 3 . Thus 
both the bias and the energy spin-splitting are of the 
same order of magnitude as the ones used in Section ITVl if 
measured in units of T. As expected the current reaches 
a steady-state in the long time limit, since the bias is 
constant and Hb ias does not have bound states^^ 9 How- 
ever, the weak links between the QD and the leads allows 
us to study a pseudo-stationary regime, in which the cur- 
rent displays well-defined quantum beats. 

At V s f = 0, the dominant frequency of the spin- 
polarized current if is = \Ul/2 + E z /2\ w 0.37 while 
1 1 oscillates with a dominant frequency cjJ = \Ul/2 — 
E z /2\ « 0.33, see panel b) of FigEl As in the WBLA, 
the difference between Wf and wj, leads to quantum beats 
in both I tot and I sp i n [see panel a) of Fig [7] and Fig[8]. 
For nonvanishing V s f the two fundamental frequencies 

lu® renormalizes as lu® — > uj a = \Ul/2 + a yjE%/4~+V£f\. 

The system Hamiltonian is no longer diagonal along the 
quantization axis z and I a acquires the second frequency 
uj—a besides the original (but renormalized) one u> a . We 
also observe that as V s f increases, the amplitude of the 
quantum beats in I to t and I sp in is suppressed, similarly 
to what happen by treating the leads in the WBLA. 

We would like to end this Section by pointing out that, 
for leads with a finite bandwidth the current might dis- 
play extra oscillation frequencies corresponding to tran- 
sitions either from or to the top/bottom of the bands, 
an effect which cannot be captured within the WBLA. 
We have investigated this scenario by changing the input 




time CO 



FIG. 7: (Colour online) I to t(t) and I to t{u) for V sf = 
[panels a)-b)], 0.05 [panels c)-d)], 0.1 [panels e)-f)] obtained 
with a bias Ul = 0.7, Ur = applied to leads of length 
N = 450. The discrete Fourier transform is calculated us- 
ing 280 equidistant points of Itot(t) — itot(oo) with t in the 
range (70,350). The rest of the parameters are V = 0.03, 
£da = 0{t)UL/2 + aE z /2, E z = 0.04 and zero temperature. 




FIG. 8: (Colour online) I a pm{t) for the same system param- 
eters of Fig[7] 



parameters in such a way that the transitions from the 
resonant level of the QD to the bottom of the band are 
energetically favored. We set Ur = 0.4, Ul = 0, E z = 0, 
E F = -1.8, e d = -0.5, V = 0.02 and V sf = 0. The 
corresponding (spin-independent) current and its Fourier 
spectrum are shown in FigjS] One can clearly see two 
well defined peaks at energies u>\ w 1.3 and 0J2 ~ L5. 
As expected, one of these frequencies corresponds to a 
transition from the resonant level to the Fermi energy of 
the L lead, i.e. \Ep — ej\ = 1.3. Transitions between the 
resonant level and the Fermi energy of the R lead, i.e. 
\Ef + Ur — Sd\ = 0.9, are strongly suppressed since the 
current is measured at the L interface. Thus the second 
peak has to be ascribed to a transition which involves 
the bottom of the band Eb = —2, specifically the tran- 
sition of energy \Eg — £d\ — 1-5. This kind of features 
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in the Fourier spectra of the transients points out the 
limitations of the WBLA and might be experimentally 
observed in QD's connected to narrow-band electrodes. 




time 

FIG. 9: Itot{t) in units of 10" 6 for U R = 0.4, U L = ap- 
plied to leads of length TV = 250. The inset shows Itot(ui) 
in units of 10 -5 calculated using 138 equidistant points of 
Itot{t) — Jtot(oo) with t in the range (12,150). The rest of the 
parameters are E z — 0, Ef = —1.8, £d = —0.5, V = 0.02, 
V 3 f — and zero temperature. 



B. Normal case: engineering the spin-polarization 

In this Section we exploit the full knowledge of the 
time-dependent response of the system in order to engi- 
neer the spin-polarization of the total current. In par- 
ticular we are interested in maintaining the polarization 
ratio 

f.y. _ o Itot(t)I S pin(t) (0A \ 

r(t) ~ iL(t) + i* pm (t) (24) 

above some given value in a sequence of time windows of 
desired duration. The quantity r is ±1 for fully polarized 
currents and zero for pure charge or spin currents. 

In FigQJJ] we show the time-dependent currents If (t) 
and I±(t) as well as Itot{t) and I sp i n (t) at zero tem- 
perature for different values of E z = 0.1, 0.6, 2 and 
U L = -U R = 1, V = 0.2, and V sf = 0. In the long-time 
limit the values of r increases steadily as E z is increased, 
and remains below 10 _1 for E z < 2. The polarization ra- 
tio can, however, be much larger in the transient regime 
and its maximum has a nontrivial dependence on E z . At 
finite values of E z there is an initial unbalance of the spin 
up and down densities. This implies that for small times 
1 1 is suppressed and delayed with respect to Jj, since 
spin j electrons can freely flow while there cannot be any 
spin I flow before the initial spin j electrons have left the 
QD»i2. This is the so called Pauli blockade phenomenon 
and is responsible for a recoil of I± during the rise of I-\ . 
In panel a)-b) of FigfTUJwe consider the case E z = 0.1. 




time time 

FIG. 10: (Colour online) Currents in units of 10 -2 for Ul — 
— Ur = 1 applied to leads of length N — 120 for different 
values of E z — 0.1 [panels a)-b)], 0.6 [panels c)-d)], 2 [panels 
e)-f)]. The rest of the parameters are V = 0.2, V 3 f = and 
zero temperature. 

Both it and I\ overshoot their steady-state value and os- 
cillate with rather sharp maxima and minima during the 
initial transient, see panel a) of FigQJjJ However, the os- 
cillation frequencies are very similar (small E z ) and the 
time at which I-\ has a maximum in correspondence to 
a minimum of I± occurs when the amplitude of their os- 
cillations is already considerably damped. On the other 
hand, for E z — 0.6, see panel c)-d) of FigQJjJ there is 
a synergy between the Pauli blockade phenomenon and 
the frequency mismatch. At small t such synergy gener- 
ates large values of the ratio r w 1 despite r(t — > oo) is 
vanishingly small. By increasing further the value of E z 
we observe a long overshoot of if while J; oscillates with 
high frequency, see panel e)-f) of FigfTUl where E z = 2. 
The ratio r(t) is a smooth decreasing function of time for 
t > 3 and approaches the value of about 10 -2 for t — > oo. 
This behavior differs substantially from the one obtained 
for E z — 0.6 where r(t) has a sharp peak for < t < 3 
and is very small otherwise. Below we show how one can 
exploit this kind of transient regimes to maintain persis- 
tently a large value of r. 

We consider the optimal case E z = 0.6 and apply a 
pulsed bias with period T and amplitude Ul = —Ur = 1 
in lead L and R respectively. The period T of Hu as is tai- 
lored to maintain the polarization ratio r(t) above ~ 0.5 
in a finite range of the period. For time-dependent biases 
Eg. (|2"2")1 has to be generalized, as the evolution operator is 
no longer the exponential of a matrix. We discretize the 
time and calculate the lesser Green's function according 
to 

G<(t n ,t n )*e- iH »^ At G<(i B _ 1 ,t n _ 1 )e w -( t ») At , 

(25) 

where t n = nAt, At is the time step, n is a positive 
integer and G<(0,0) = if(H ). 

In Fig llll we plot two time-dependent responses for 
T = 6 and T = 8. The pulsed bias produces an alternate 



9 



Itot and I sp in whose amplitude depends on T. We note 
that the amplitude of I tot is of the same order of mag- 
nitude of the steady state value hotit — ► oo) attained in 
FigHSl panel d) for constant bias. On the contrary the 
amplitude of I sp i n is two orders of magnitude larger than 
the corresponding steady-state value. The polarization r 
cannot be maintained as large as r w 1 (maximum value 
of r during the transient) due to an unavoidable damp- 
ing. However, the value of r is above 0.5 in a time window 
of 1.4 for T = 6 (with a maximum value r ~ 0.75) and 
2.6 for T = 8 (with a maximum value r ~ 0.9) in each 
period, see FigfTTI 




time time 

FIG. 11: (Colour online) Itot{t) and I sp i n {t) in units of 10 -2 
for a pulsed bias of period T = 6 (top-left panel) and T = 8 
(bottom-left panel) of amplitude Ul = — Ur — 1 applied 
to leads of length N = 120. The polarization ratio r(t) is 
displayed in the right panel. The rest of the parameters are 
the same as in Fig llOl The dotted lines in the left panel 
denote the pulsed bias in the left lead in units of 10 -1 . The 
numerical calculations have been performed with At = 0.1. 



C. Ferromagnetic case: TMR 

The spin-dependent band-structure of the leads in- 
troduces another frequency dependence in the currents. 
Both 7j and I± display a coherent oscillation of frequency 

^1,2 = \hi - h 2 \ = ' El/A + Vfj, where h\,2 arc the 
eigenvalues of the isolated QD. Such behavior cannot be 
observed in the normal the model is diagonal 

along the quantization axis £ of Eq. (|16[) and no tran- 
sitions between states of opposite polarization along £ 
can occur. In Fig |12l we display the spin-up current at 
zero temperature in the P configuration for Ae = 1.5, 
Ul = 1, U R = 0, e d<J = 6{t)U L /2 + aE z /2, E z = 0.04, 
V = 0.02, and V s f = 0.1. Besides the frequencies lol\ = 
\E F + U L -hi\n 0.6 and uj L2 = \ E F + U L - h 2 \ « 0.4 
already discussed in Section IIII Al one can see the ap- 
pearance of the new frequency io\ i2 ~ 0.2. 

Next, we study the steady-state regime in both the P 
and AP configurations and calculate the TMR for differ- 
ent values of the bias voltage U, the band spin splitting 
Ae and the spin- flip energy V s f. Analogies and differ- 
ences with the case of wide-band leads will be discussed. 

In Fig[T3] we display the contour plot of the TMR 
in the parameter space spanned by the bias voltage 
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time 

FIG. 12: J T (t) in units of 10 -4 for U L = 1, Ur = ap- 
plied to leads of length N = 350. The inset shows the dis- 
crete Fourier transform I^(u}) calculated using 265 equidistant 
points of Ii(t) — Jf (oo) with t in the range (35,300). The rest 
of the parameters are e dtT = 6(t)Ui,/2 + aE z /2, E z = 0.04, 
V — 0.02, V a f =0.1 and zero temperature. 



U = Ul = —Ur and the band-spin-splitting Ae at in- 
verse temperature (3 = 100, V = 0.5, e^ = and for 
V s f = 0, 0.25, 0.5. In the left panel we show the TMR 
for V s f — 0. We can see that despite the dot-leads link 
is symmetric, there is a finite region at small bias and 
intermediate Ae in which the TMR < 0, although rather 
small (« —0.05). This is a new scenario for the TMR in- 
version and stems from the finite bandwidth of the leads 
(we recall that TMR > for V sf = in the WBLA). The 
largest positive value of the TMR (TMR pa 0.7) occurs 
for large magnetization and bias, as expected. 

As V s f is increased (central and right panels of Fig [T5|) 
the region of positive TMR widens, which is an opposite 
behavior to the one in the WBLA. However, we note that 
the largest value of negative TMR occur at small U and 
large Ae (TMR pa -0.2, see right panel of FigIT3]) and 
that the TMR changes sign abruptly as U is increased, 
a feature which is in common with the WBLA. The pos- 
itive values of the TMR reduce with respect to the case 
V s f = 0. This property is due to the presence of spin- 
flip scatterings which close conducting channels in the P 
configuration and open new ones in the AP configuration, 
thus suppressing the difference Ip — Iap ■ 

Finally we have investigated the dependence of the 
above scenario on temperature. It is observed that the 
qualitative picture described in Fig |13l survives down to 
(3 pa 10. Increasing further the temperature the TMR 
behaves similarly to the WBLA. 



V. SUMMARY AND CONCLUSIONS 

The ultimate goal of future QD-based devices is the 
possibility to generate spin-polarized currents, control 
their spin-coherence time, and achieve high TMR after 
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FIG. 13: (Colour online) Contour plot of the TMR at the 
steady-state in units of 10~ 2 as a function of U and Ae for 
different values of V sf = (left), 0.25 (center), 0.5 (right). 
The boundary TMR = is displayed with a white line. The 
remaining parameters are V = 0.5, Eda = 0, (3 = 100 and the 
length of the leads is N = 100. 

the application of high-frequency signals. This calls for 
a deep understanding of the time-dependent responses in 
these systems. 

In this paper we have calculated spin-dependent out-of- 
equilibrium properties of lead-QD-lead junctions. Realis- 
tic transient responses are obtained within the partition- 
free approach. The time-dependent current is calculated 
for QDs connected to ferromagnetic leads and in the pres- 
ence of an intra-dot spin flip interaction. This requires 
the propagation of a two-component spinor. 

For ID leads, we evolve exactly a system with a finite 
number N of sites in each lead. If N is sufficiently large, 
reliable time-evolutions are obtained during a time much 
larger than all the characteristic time-scales of the infi- 
nite systemj 25 i 26 i 27 By comparing our results against the 
ones obtained with leads of infinite length t 24 i 28 we have 
verified that our method is accurate and robust, beside 
being very easy to implement. 

We have solved analytically the time-dependent prob- 
lem in the WBLA and derived a closed formula for the 
spin-polarized current. Such formula generalizes the one 
obtained in the spin-diagonal case and has a transparent 
interpretation! 10 i 30 We stress, however, that within the 
WBLA, transitions involving the top or the bottom of the 
leads band are not accounted for. The latter may be rel- 



evant to characterize, e.g, the coherent beat oscillations 
when the device is attached to narrow-band electrodes. 

Furthermore we have shown how to engineer the 
transient response of the system to enhance the spin- 
polarization of the current through the QD. This is 
achieved by controlling parameters like, e.g., the external 
magnetic field, the transparency of the contacts, and im- 
posing a pulsed bias of optimal period. It is shown that 
by exploiting the synergy between the Pauli- blockade 
phenomenon and the resonant-continuum frequency mis- 
match, one can achieve an AC spin-polarization two or- 
ders of magnitude larger than the DC one. 

We also have employed the Stoner model to de- 
scribe ferromagnetic leads and computed the steady-state 
TMR. We have found a novel regime of negative TMR, 
in which the geometry of the tunnel junction is not re- 
quired to be asymmetric and a finite intra-dot spin-flip 
interaction turns out to be crucial. For any given V s f 
there is a critical value of the ferromagnetic polarization 
above which the TMR is negative. The magnitude of the 
TMR is very sensitive to temperature variations and the 
TMR inversion phenomenon disappears as [3 approaches 
the damping-time of the system. 

We would like to stress that our approach is not limited 
to ID electrodes and can be readily generalized to inves- 
tigate multi-terminal devices consisting of several multi- 
level QDs. Finally, owning to the fact that the prop- 
agation algorithm is based on a one-particle scheme, it 
prompts us to include electron-electron interactions at 
any mean-field level or within time-dependent density 
functional theory^ 
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